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Abstract 

This paper is dedicated to the study of viscous compressible barotropic fluids in 
dimension N = 2. We address the question of the global existence of strong solutions 
with large initial data for compressible Navier-Stokes system and Korteweg system. 
In the first case we are interested by slightly extending a famous result due to V. 
A. Vaigant and A. V. Kazhikhov in [32] concerning the existence of global strong 
solution in dimension two for a suitable choice of viscosity coefficient — fi > 

and X(p) = XpP with f3 > 3) in the torus. We are going to weaken the condition 
on (3 by assuming only f3 > 2 essentially by taking profit of commutator estimates 
introduced by Coifman et al in [5] and using a notion of effective velocity as in |32j . 
In the second case we study the existence of global strong solution with large initial 
data in the sense of the scaling of the equations for Korteweg system with degenerate 
viscosity coefficient and with friction term. 

1 Introduction 

The motion of a general barotropic compressible fluid with capillary tensor is described 
by the following system, which can be derived from a Cahn-Hilliard free energy (see the 
pioneering work by J.- E. Dunn and J. Serrin in [8] and also in [21 El 111]): 

\ d t p + d\v{pu) = 0, 

[ d t (pu) + d'w(pu g) u) - div(2p(p)Du) - V(X(p)divu) + VP(p) + apu = divK, 

(1-1) 

where divK is the capillary tensor which reads as follows: 

divK = V( P K(p)Ap + + P k (p))|V p| 2 ) - div(K(p)Vp® Vp). (1.2) 

The term divK allows to describe the variation of density at the interfaces between two 
phases, generally a mixture liquid-vapor. P is a general increasing pressure term that we 
assume in the sequel under the form P(p) = bp' 7 with b > and j > 1, apu is a friction 
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term with a > (see [25]). D(u) = ^[Vit +* Vu] being the stress tensor, p and A are the 
two Lame viscosity coefficients depending on the density p and satisfying: 

p > and 2p + NX > 0. 

In the present paper, we are interested in dealing with two different situations: 

• The case of the compressible Navier-Stokes system where we assume no capillarity, 
n(p) = and where p,(p) = 1 is a constant and X(p) = Xp@ with (3 > 2. 

• The case of Korteweg system with the viscosity coefficients and the capillarity 
coefficient such that: 

pip) = pp X(p) = and nip) = — with k > 0, p > 0. 

p 

In the first case we would like to extend the famous result of global strong solution 
in two dimension in the torus discovered by V. A. Vaigant and A. V. Kazhikhov in 
[52"] . Indeed in [32], the authors assume that X(p) = Xp@ with j3 > 3 and A > and 
p(p) = 1. Let us emphasize that a such choice on the viscosity coefficients allows to 
exhibits two different phenomena; the first one concerns the notion of effective velocity 
introduced by Lions in [23] which is crucial in this context for getting a priori estimates 
on the divergence and the rotational of the velocity; the second reason to choose such 
coefficient concerns the possibility to obtain L^?(L P (T 2 )) estimates for any p > 1 on the 
density and any T > 0. Indeed the viscosity coefficient X(p) offers enough weight in order 
to derivate such estimates on the density (see the pi 115 "Second a priori estimate for 
the density" in [32]). In the first part of this paper we wish to improve this result by 
assuming only X(p) = Xp 13 with (3 > 2 and A > 0. The key point will consist in using 
commutator estimates for dealing with term of the form [Rij, Uj](pui) (we refer to |6j for 
such estimates but also to ?? where Lions prove global existence of weak solution for 
compressible Navier-Stokes equations by introducing this kind of ingredient). 
In the second case we are interested in proving global existence of strong solution for the 
Korteweg system with friction term when the physical coefficients verify: 

k(p) = —, k = p 2 and b = ap, (1.3) 

with p > 0. This system without friction has been widely studied this last year in 
particular concerning the existence of global weak solution and global strong solution with 
small initial data. We refer in particular to the following works [^[71 115 1 [16 ] 117 1 [18l 1191 120j. 
Let us start with explaining a other notion of effective velocity used in particular in [15] 
which allows us to simplify the system (|1.6p . Indeed by computation (see [15] ) , we obtain 
the simplified system: 

d t p + div(pv) - pAp = 0, ^ 
pdtv + pu ■ Vv — div(/ip Vv) + a p v = 0, 

with v = u + pV In p the effective velocity. For more details on the computation, we 
refer to |16] , When we write the system (|1 .4[) in function of the momentum m = pv, the 
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system reads as follows: 

dtp + divm — pAp = 0, 

m (1.5) 
dtvn + div( — ®m)-/iAm + rm = 0, 

P 

In particular we observe that (p, — uVlnp) when: 

d t p - pAp = 0, 
is a particular global solution of the system (|1.6p . 

Remark 1 Let us mention that we can choose initial density which admits vacuum that 
we wish. In general it is not always possible to obtain global strong solution with initial 
density close from the vacuum. 

In the sequel we will be interested in working around this global particular solution (see 
remark [1]) in order to prove the existence of global strong solution with large initial 
data on the irrotational part. More precisely we shall wish obtain global strong solution 
in critical space for the scaling of the equations. Let us briefly recall the notion of 
invariance by scaling of the equation and by what we mean by supercritical smallness 
on the data. By critical, we mean that we want to solve the system (jl.6p in functional 
spaces with invariant norm by the natural changes of scales which leave (jl.6p invariant. 
More precisely in our case, the following transformation: 

(p(t, x),u{t, x)) — > {p{l 2 t, lx),lu(l 2 t, lx)), I G R, (1.6) 

verify this property, provided that the pressure term has been changed accordingly. In 
particular we can observe that H 2 x H 2 is a space invariant for the scaling of the 
equation, more generally such Besov spaces: 

JV N y 

( Po -p)eB; A , uo€B£ tl , 

with (p,pi) G [1, +oo[ are also available. 
1.1 Results 

Let us state the two main result of this paper. The first one is an improvement of the 
results of Vaigant and Kazhikhov [32] . 

Theorem 1.1 Let us assume the following hypothesis on the viscosity coefficients: 

p(p) = 1 and \{p) = Xp 13 with /3 > 2. 
Let n e H 2 (T 2 ), p e W^ q {T 2 ) with q>2 and: 

< c < po(x) < m < +00 Vx G T 2 , 
then it exists a unique global strong solution to h 1.0(1 such that: 

u G W 2 '\Q T ) and p G W^(Q T ) VT > 0, 
with Q T = (0,T) x T 2 . 
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Remark 2 As mentioned above, the main point compared with the result of \32$ concerns 
the fact that we can improve the range of j3 by assuming only (3 > 2. It would be possible 
also to improve the regularity condition on the initial data by working in Besov space 
invariant for the scaling of the system, but it is not the object of this paper. 

We are going to give your second result on Korteweg system with supercritical smallness 
condition on the initial data; before let us give the following definition: 

Definition 1.1 We set q = p — 1, m = pu and p = h. 

In the following we are dealing with the euclidian space with N > 2. Let us give 
our main result on the Korteweg system where we prove the existence of global strong 
solution with large initial data on the irrotational part. 

Theorem 1.2 Suppose that we are under the conditions Assume that uiq G B 9 2 1 

JV 

and qo G -B 2 2 i with ho > c > 0. Then there exists a constant sq depending on ^ such 
that if: 

\\mo\\ 4-1 < £ o, 

then there exists a unique global solution (q, m) for system hl.5\) with h bounded away 
from zero and, 

h € C(R + , Bfa) n L 1 (M + , Bfc ) and m G C(R + ; B£ x ) n L 1 (M+, B 2 2 A n B 2 2 ;1 ). 

Let us give the plane of this paper, we shall remind in section [2] some auxiliary results of 
Gagliardo-Nirenberg's inequality and in section Othe Litllewood-Paley theory. In section 
H] and section O we will prove different a priori estimates on the density and the velocity 
which show the theorem II. H We will conclude in section [U] by the proof of theorem II. 21 

Notation 

In all the paper, C will stand for a harmless constant, and we will sometimes use the 
notation A < B equivalently to A < CB. 

2 Auxiliary Assertions 

We are going to recall some lemma which are also present in [32J and that we prefer to 
state for the sake of the completnes. 

Lemma 1 Let f2 G M. N be an arbitrary bounded domain satisfying the cone condition. 
Then the following inequality is valid for every function u G W 1,m {^l), J n udx = 

IMU«(n) ^ C i\\^ u \\lrn (n) \\u\\ l -* nv (2.7) 
where a = , r , q i , moreover if m < n then q G \r, mn 1 for r < mn and q G [ m " , rl 

r m n 

for r > . If m > n then q G [r, +oo0 is arbitrary; moreover if m > n then equality 
{2. 7\) is also valid for q = +oo. 
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Inequality (|2.7p is a particular case of the more general inequalities proven in [2TJ [TO] . 
Let us mention that an inequality of the form (|2.7p is valid for the function of class 
W 1,m (Q) when M = -py Jq u dx is not null. It suffices to consider v = u — M and apply 
inequality (j2.7|) to the function v. We obtain then the inequality 

||«IU»(n) < C2(||Vu||2 m(n) ||u||^ ( Q c) + ||n|| z i ( n)), (2.8) 



Lemma 2 Let £1 £ M 2 be an arbitrary bounded domain satisfying the cone condition. 
Then every function u G W 1,m (£l) with J^udx = satisfies the inequality 

\\u\\ L ^. < C 3 (2 -my^\\Vu\\ L m {n)l l<m<2, (2.9) 

where C3 is a constant independent of m and the function u. 

For a proof of this inequality see |33} 130] . The exact constant in inequality (|2.9p is 
obtained in the article [30] . 

Lemma 3 Let S] £ I 2 k an arbitrary bounded domain satisfying the cone condition. 
Then for an arbitrary number e, 1 > 2e > 0, every function h G W ' m + s (Sl), m > 2, 
1 > J > 0, satisfies the inequality 

\\h\\v*» {n) < C4(Nki(fi)+ mi HV^fX, l^ll^(i-.)(n)). ( 2 - 10 ) 

L m +* (Q) v ; 

where s = (1 — e) m _^f_ £ ) C4 is a positive constant independent of m,e,8 and the 
function h. 



3 Littlewood-Paley theory and Besov spaces 

Throughout the paper, C stands for a constant whose exact meaning depends on the 
context. The notation A < B means that A < CB. For all Banach space X, we 
denote by C([0, T],X) the set of continuous functions on [0, T] with values in X. For 
p € [1, +00], the notation L p (0, T, X) or Lj,(X) stands for the set of measurable functions 
on (0,T) with values in X such that t — > \\f(t)\\x belongs to L p (0,T). Littlewood-Paley 
decomposition corresponds to a dyadic decomposition of the space in Fourier variables. 
We can use for instance any <p G C°°(M N ), supported in C = {£ G M^/f < |£| < |} such 
that: 

5>(2- J o = i if e/o. 

Denoting h = J 7 " 1 ^, we then define the dyadic blocks by: 

A lU = (p{2' l D)u = 2 lN / h(2 l y)u{x - y)dy and S t u = V A k u . 

Formally, one can write that: 

u = ^ A fc u. 

This decomposition is called homogeneous Littlewood-Paley decomposition. Let us ob- 
serve that the above formal equality does not hold in S (R ) for two reasons: 
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1. The right hand-side does not necessarily converge in <S 



!>N\ 



2. Even if it does, the equality is not always true in S'(M N ) (consider the case of the 
polynomials). 

3.1 Homogeneous Besov spaces and first properties 

Definition 3.2 For s£l, p£ [1, +oo], q G [1, +oo], and u G S' (R N ) we set: 

The Besov space B^ q is the set of temperate distribution u such that ||u||b<> < +oo. 

Remark 3 The above definition is a natural generalization of the nonhomogeneous Sobolev 
and Holder spaces: one can show that B^^ is the nonhomogeneous Holder space C s and 
that B2 2 is the nonhomogeneous space H s . 

Proposition 3.1 The following properties holds: 

1. there exists a constant universal C such that: 
C -1 ||u||rs < ||Vu|| R s-i < C||u||rs . 

11 11 p,r — 11 "-Op,r — n-^PjT 

2. If Pl < P2 and n < r 2 then B s puri B^ 1/pi ~ 1/pa) . 

' 1 1 

3. Bp ri Bp r if s > s or if s = s and n < r. 

Let now recall a few product laws in Besov spaces coming directly from the paradiffer- 
ential calculus of J-M. Bony (see [?]) and rewrite on a generalized form in [?]. 

Proposition 3.2 We have the following laws of product: 

• For all s£l, (p, r) G [l,+oo] 2 we have: 

\\uv\\ B ^ r < C(||«||L°°lb||B| ir + |M|L°o||u|| B ^ r ) ■ (3.11) 

• Let (p,pi,p 2 ,r,\i,\ 2 ) G [l,+oo] 2 such that:± < ^ + ^, p\ < \ 2 , p 2 < M, \ < 
p\ + and ^ < ^ + j^. We have then the following inequalities: 

if si + 82 + iVinf(0, l-±-X) >Q)8l+ N < K anda2+ N < K then: 

\\uv\\ _ n( j_ + j__^ < \W\ B ; lr \\v\\ B;iao , (3.12) 

when Sl + ^ = £ (resps 2 + ^ = f J we replace ||«|| B ;i r IMIb^ ( res P \\ v \\b%,J 
KU Mb ^ (re$P II^H^.oonLoo;, if 8! + g = f and s 2 + £ = g 



by ||«|| B ;i 1 IMIb^ ^ ||w||^a >00 nLooA */ «i + £ = f and 82 + £ = f we taA;e 
r = 1. 

// Sl + S2 = ; Sl G(f -g] amli + i < 1 then: 

-£*p,oo 
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If \s\ < — for p > 2 and — ^ < s < — else, we have: 
p p p ' 



N 



N 



\uv\\bs , < C||u]|bs llull jv 



(3.14) 



Remark 4 In the sequel p will be either p\ or p2 and in this case j = ^- — ^ i/pi < P2, 

res p\ = h~p\ tfpi^pi- 

Corollary 1 Let r £ [1, +oo], 1 < p < pi < +oo and s suc/i that: 
v Pi ' Pi I J P Pi — ' 

. s g ( — — + N(- + — — 1), — ) if - + — > \, 

V pi Vp pi <" pi / J p pi ' 

JV 

then we have if ' u £ Bp and v £ B^ t00 n L°° : 

\\uv\\b" < C||«||rs 1 1 17 1 



The study of non stationary PDE's requires space of type L p (0,T, X) for appropriate 
Banach spaces X. In our case, we expect X to be a Besov space, so that it is natural to 
localize the equation through Littlewood-Payley decomposition. But, in doing so, we ob- 
tain bounds in spaces which are not type L p (0, T, X) (except if r = p). We are now going 
to define the spaces of Chemin-Lerner in which we will work, which are a refinement of the 
spaces Lj,(Bp ). 



Definition 3.3 Let p £ [l,+oo], T G [l,+oo] and si G R. We set: 

i«iiz^ ( B-) = (E 2 ' rsi ii A ^)ii^(^))"- 



We then define the space Lj,(Bp^ r ) as the set of temperate distribution u over (0,T) x R N 
such that \\u\\yp rFl si \ < +oo. 

We set Ct{B*\) = L^(Bp l r ) t~) C([0,T], B^ r ). Let us emphasize that, according to 
Minkowski inequality, we have: 

IMIz^Km ^ IMI;4« r ) if r ^ p> W u h? T {B% r ) ^ ll«ll^« r ) if r < p- 



Remark 5 It is easy to generalize proposition VS. "A to L P r ,(Bp 1 r ) spaces. The indices s\, 
p, r behave just as in the stationary case whereas the time exponent p behaves according 
to Holder inequality. 

In the sequel we will need of composition lemma in Lj,(Bp r ) spaces. 
Lemma 4 Let s > 0, (p,r) £ [l,+oo] and u £ L p T [B s pr ) n Lf{L°°). 
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1. Let F G W^ ] C +2,C °(R N ) such that F(0) = 0. Then F(u) G L p T (B s pr ). More precisely 
there exists a function C depending only on s, p, r, N and F such that: 



W F ( u )h^ r ) ^ c (\\ u h¥(L°°))\\ u k? T (Bi r y 

2. Ifv,u€ L p T {B s pr ) n Lf{L°°) and G G W^ +3 '°°(R N ) then G{u) - G(v) belongs to 
L^(Bp) and there exists a constant C depending only of s,p,N andG such that: 

\\G(u) - G(v)\\^ {B s r) < CflMlioo^oo), - uh^B'^i 1 + \\ u \\l™(l°°) 

+ IMU§?(L°°)) + \\V - u\\L ¥ (L<x>)(\\u\\z? T (BS )r ) + IMIZ£(B3 ir )))- 

Now we give some result on the behavior of the Besov spaces via some pseudodifferential 
operator (see [?]). 



Definition 3.4 Let m G M. A smooth function function f : M — > M is said to be a S r 
multiplier if for all muti-index a, there exists a constant C a such that: 

V(eR N , \d a f(0\<C a (l + \m m ' H - 



Proposition 3.3 Let m G R and f be a S m multiplier. Then for all s £ M and 1 < 

p,r < +oo the operator f(D) is continuous from B pr to Bp~ m . 

Let us now give some estimates for the heat equation: 

Proposition 3.4 Let s G R, (p, r) G [l,+oo] 2 and 1 < p2 < p\ < +oo. Assume that 
uq G r and / G L P ^{Bp^ r 2+2 ^ P2 ). Let u be a solution of: 

j d t u - /tAu = / 
[ u t =o = u . 

Then there exists C > depending only on N, p,, p\ and p2 such that: 

IHIz«(s;+ r 2 ^) ^ c ^M\b 1)T + ^' l \\f\\~ L P2 {K:r 2 + 2 /P2) ) . 

If in addition r is finite then u belongs to C([0, T], Bp r ). 

4 Proof of theorem 11.11 

In the sequel we shall work on the torus O = T 2 . Let us start with recalling the energy 
estimate, when we multiply the momentum equation we get: 

(p\u\ 2 (t, x) + 7r(p)(i, x))dx + / / \Vu\ 2 (s,x)dsdx 
n Jo Jn 

+ 11 (1 + \{p)(s,x))(dwu) 2 {s,x)dsdx < / {p {x)\u (x)\ 2 + U(p )(x))dx 
Jo Jn Jn 



(4.15) 
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with 7r defined as follows: 

k(p) = a {—^—- r(P 7 - P) ~ P + !) for 7 > 1. 
7-1 

Let us recall that P' (p) = pir" (p) what implies by convexity that 7r(p) > 0. Finally as 
we assume that: 

c i= [ {t;Po( x )\ u o( x )\ 2 + 7T (po)(x) + Po(x))dx 
Jn 1 

is finite, we obtain at least formally (if p and u are enough regular for performing inte- 
gration by parts) by energy estimate (|4.15p and via the transport equation that: 

(p\u\ 2 (t, x) + ir(p)(t, x) + p(t, x))dx + I I | Vu| 2 (s, x)dsdx 
v. Jo Jn 

t 



+ I / (1 + \(p)(s,x))(divu) 2 (s,x)dsdx < C 



(4.16) 

l 

io Jn 

Let us now explain how to get £ 2 ((0, T) x fi) estimates on u, we are going to follow Lions 
in [23] p4. Indeed by the momentum equation we have: 



pu(t,x)dx\ = | / p u (x)dx\ < \\poUo\\ L i(Q). 
n Jn 

Next we use the Poincare-Wirtinger inequality and we have: 



p(t,x)[u(t,x) - / u(t,y)dy]dx\ <C\\p(t,-)\\ L 7\\Vu\\ L 2( Q y 
n Jn 

Hence for all t > 0: 

u(t,x)dx\ < 1 (\\poUo\\Li(n) + C\\p(t,-)\\ L -f\\Vu\\ L 2 {n) ). 
Un Po ax ) 

We conclude by Poincare-Wirtinger inequality which implies that | u(t, x)dx\ + \\Vu\\i i 2^ 
is an equivalent norm to the usual one in H 1 ^). □ 

Now we are just going to explain where in the proof we can slightly improve the range 
of the coefficient j3 in |32| . One of the main point of the proof in [32] consists in get- 
ting a priori estimates on the density in L°°(L P (T 2 )) for any p > 1. This is possible 
due to the viscosity coefficient X(p) = p 13 which provide such estimate at least if j3 is 
large enough. Let us follow the arguments of the proof of [32] and explain where by 
commutators estimates we can weaken the hypothesis /3 > 3. 

5 A priori estimates on the density and the velocity 



First as in [32], we are going to recall some estimates for solutions to the following two 
Neumann problems: 

A£ = div(pu), / £dx = 0, d Xl t,\ Xl=0 , Xl=1 = d X2 £\ X2=0tX2= i = 0. (5.17) 
Jn 
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At] = div(div(pu ®u)), / rjdx = 0, , d xl r)\ xl =o !X1 =i = d X2 rj\ X2=0 ^ X2=1 = 0. (5.18) 

Jn 

Therefore by [23] we have solution to the problems (|5.17p and (j5. 18|) . whereas the esti- 
mates for singular integrals in [28] provide the following inequalities: 

||V(A) _1 div(pti)|| L 2m < m||/)ti|| L 2m, 1 < m < +oo, 
||V(A)- 1 div(pu)|| L 2-r < \\pu\\ L 2- r , 1 > 2r > 0, 
\\Ri,j (puiUj)W L 2m ^5 m\\pu ® u\\ L 2m, 1 < m < +oo, 

Here we have roughly written £ = (A) _1 div(pw) and rj = Rij(puiUj) (with the summation 
notation) . 

By Holder's inequalities we obtain: 

||V(A) _1 div(pu)|| i 2m < m||p|| 2mk \\u\\ L 2mk , 

L k- 1 

HvCAj^divMII^- £ llp|| L a-|lvHli, ( 5 - 19 ) 

\\Rij(pUiUj)\\ L 2 m <m\\p\\ 2rnk\\u\\ 2 L4mk , 

where k > 1, m > 1 and r > 1, 1 > 2r > 0. 
From the estimate of lemma [T]l3l we obtain: 



IMIl 2 '™ ^ m2 II^^IIl 2 ; m> 2, 

||(A) _1 div(pn)|| L 2m < mi || V(A) _1 div(pit) || _2m_ , m > 2, 



We set now: 



(5.20) 



= / (curhi 2 (i,x) + (2 + A(p))divu 2 (i,x))cfe, 

we obtain then by using (|5,20p , (|5.19p with r = ^q-j- and the energy inequality (|4.16p : 

\\{AY x ^{ P u)\\^ m <wh\\p^ Lm , m>2, (5.21) 

Similarly we have: 

||V(A)" 1 div(pu)|| i 2m < m%k^(ip(t))5 \\p\\ m, m > 2, k > 1, 

Lfc " (5.22) 
||^j(/9n^Uj)||^2m < m &</?(f)||p]| 2mfe , m > 2, fe > 1. 

5.1 Gain of integrability for the density 

Following |32] the plan of the proof of [32], we are interested in getting a gain of inte- 
grability on the density. We follow here the method of Lions in |24] to get a gain of 
integrability on the pressure and the argument developed in [32J. Apply the operator 
(A) _1 div to the momentum equation, we obtain: 

d 

— (A)- x div(H + [Ri^u^pui) - (2 + \(p))dwu 

(5.23) 



+ P{P) ~ j4 J( p (p)^ x ) ~ ( 2 + Kp))^u)dx = 0. 
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We will set in the sequel: 

B = (2 + A(p))divu-P(p). (5.24) 
Next if we renormalize the mass equation we have: 

d t 6(p) + u ■ V9(p) + p6'(p)divu = 0. 

where we have set: 

0(p) 



[ P -(2 + X(s))ds = 21np + \{pP - 1). 
Ji s P 



(5.25) 



Finally we get the following transport equation: 

J^[(A)-Miv(H + 0(p)] + u-V[(A)- 1 dw(pu) + e(p)] + [R ij ,u j }(pu i ) 

+ Pip) ~W\J n t P(p)( *' x) " (2 + x (p)) d[vu \ dx = °. 

Denote by / the function: 

/(t,x) = max(0, (A) _1 div(/9it) + 9{p)) 

and multiply the equation (|5.25p by the function pf 2m ~ l with m £ N and m > 4 and 
integrate over , we obtain: 

^~T+ I Pf mdx + [ pP(p)f m - 1 dx+ [ [R lJ ,u 1 ](pu t )pf 2m - 1 dx 

2mdtJ n J n Jn (5 26) 

+ Bdx pf m - l dx = 0. 
Jn Jn 

As in [32] we set: 

Z(t) = ( [ pf 2m (t,x))^ (5.27) 
Jn 

Using Holder's inequality + ^ = 1, + 2m(2 ^ +1) = ^), we begin with 

estimating the term | f n [Rij,Uj](pui)pj ' 2m ~ 1 dx in (|5.26|) as follows: 



[Rij^jjipu^pf^dxl < / WR^u^pui^p^f^dx 

Jn 

< || WR^u^pu^pMh^Z 2 ™- 1 ^) (5.28) 
p||^ +1(Q) ll[%,%-](p^)ll i2ro+ ^ (fi) ^ 2m_1 W 



< 



Next we recall a result of R. Coifman, P.-L. Lions, Y. Meyer and S. Semmes in [6], which 
says that the following map: 

W hri (T N ) N x L V2 (T N ) N VF 1 ' r3 (T Ar ) Ar 
(a, 6) -)■ [a^i?^]^ 
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is continuous for any N > 2 as soon as — = — + — . Hence we obtain that [Rij, Uj\(pUi) 
belongs in (where 1 = \ + 55^55 + with fc > 1 and p = 2 - < 2) 

with the following inequality: 

\\[Rij,V,j](pUi)\\ W l,p(p) < C||V«||ia(n)||«|| 2(m+^)k \\p\\ 2(m+^)k , 

(5.29) 

< C'||V«|| 2 ^(n)||«||^a(m+^)fc (n Jp||L*»fl+i(n)i 

where we have choose k such that ^ k-i^ = ^ m fi + 1, let A; = 2m(™i)+\ — ~' ^ e ver ^ es 
that I = ~ — 5 = 2m ^j_ . Next by using lemma [2] and f|5.29j> we get: 







||[i?ij,%-](w)(t,-)|| i2m+ i < (m^\\Vu(t,-)\\ L 2 {n) \\u(t,-)\\^ m+ ^ )k \\p{t, Olli^+^n) 

+ | / [i2ij,Uj](pMi)(t,x)cix|). 

We can easily bound the last term on the right hand side by using the continuity of the 
Riez transform in L p (cv) with 1 < p < +00: 



[R ij ,u j ](pu i )(t,x)dx\ < ||p(V)lll/v(fi)IKV)llffi(n) 

! 

By (|5.2U|) and the previous inequalites we obtain finally: 



\\[Rij,Uj}(pui)(t, -)\\ 2m+ i < Cm\\Vu\\ 2 L2(n) \\p\\ L 2 m p + i {Q) + \\p(t,-)\\mn)\\u{t,-)\\ 2 H i(n)- 



(5.30) 

We have then from ([5728]) and (pTBUj) : 

[R lv u ] ]{pu l )pf 2m - 1 dx\ < (m\\p\^ +1{Q) tp{t) + Utr^^Ht,.)^^ 2 ™- 1 ^). 
Next as in [32] we get: 



Bdx I P f 2m - l dx\ < Z 2m - l {t)\\p\\l™ I ((2 + A)|divu| + P)dx 
Jn Ja 



<Z 2m ~\t){l + ^{t))\{ I (2 + X(p))dx)i) 



< z 2m -\t){\ + p(t)* + ||p||f 2ro/3+1(n) ^)5). 



Collecting all the above inequalities, we obtain: 
ft r t 



Z{t)<l + \\p\\L^n)(r)\\u\\ 2 HHn) (r)dr+ J mv(r)||p||^ +1(n) (r)dr 



(5.31) 

* 1 I 
+ / v{-r) 2 \\p\\l2mp+i {n) (T)dT. 
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As we have seen that u belongs in L 2 ((0,t), H 1 (0,)) we have: 



p\\L-<(n)(r)\\u\\ 2 HHn) (T)dT < 1. 







We obtain then: 



Next following [32] we introduce the measurable sets: 

n 1 (t) = {xen/p>2m'} and Q 2 (i) = {x € 1 (t)/6'(p) + (A) -1 div(pu) > 0}. 
We then have: 



IMI? 2 ^ + i (n) < ( ^ t p 2m/3+1 dx)^+r + lf (5.33) 
Moreover by the definition of the function 8(p), we have: 

( I p ^+i dx ^2d+i < ( f pOip^dx)^ 1 (5.34) 

Using the fact that on fii(i) \ ft 2 (i) we have < 9(p) < |(A)- 1 div(pu)|, we derive the 
following estimate: 

P 6{p) 2m dx = [ p{9{p) + (A)- 1 div(pu) - (A) _1 div(pu)) 2m (ix 



7ni(t)\n 2 (t) 



<2 2m - 1 ( pf{ P ) 2m dx+ p|(A)-Miv(pu)| 2m cfo) 
Jn 2 (t) Jn 2 (t) 



< 



+ / p\{/±)- l &\-v{pu)\ 2m dx 

2 2m (Z 2m (t)+ [ p\{A)- l d\w{pu)\ 2m dx) 
Jn 

From estimates ()5.33p and (|5.34p we deduce: 

\\pf L 2 m , + i {n) < C((m + + ( / p\{&y l d\v(pu)\ 2m dx)^f^) 

In view of estimate (|5,2ip . (|5.22j) . we have: 

H(A)~Miv(HI 2m ^< IIpIIl-^^IKA)-^^^)!! 2 ^! , 

< ||p|| w+ i (n) ((m + ^)*\\p\\ ] m+ 1 _) 2m 



< C m m m ||p l|m+1 
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20' " rn L m+ P(n) ; 



Finaly: 

2 B ()(m + l) 

l|p||? W i (n) < {Z(t)W(t) +m L 2\\ P \\l^ 1{Q) ). (5.35) 
By Young's inequality (with q = 2 ^^ 1 and p = ^Jri + m(2 ^_ 1) ), we obtain: 

^ + i (n) < (2(t) + ^m^ + ^TT + £ ||HI? W(n) )- 
By bootstrap, we get: 



\\pf L2mP+1{n) <Z(t) + m^. (5.36) 
Therefore (|5,3ip and ()5.36p give the following inequality: 

\\p\f L 2 m n + i {a) < i m ^ =l+ j Q m ^ T )\\p\\ l L^{^ T ) dT + f Q ^ T ^\\p\\l2mf J+HQ) (r)dT) 
Next by Young's inequality, we have: 

l|p||Jam/M-i ( n) £ (l+m^r+^ m<p(T)\\p\\]^ +im (T)dr + J \\pf L2m , +1{Q) (T)dr) 
Using the fact that ip(t) E L 1 (0,T) and applying Gronwall's inequality, wehave that: 

l|p||^+i (n) < C{l + mW^ +^ rop(r)|MI^+i (n) (r)dr,) 
where C depends on t. Denote: 

V(t) =m'^\\p\\ L2m/3+1(c) , t€[0,T\. (5.37) 



Then: 



i/^mf- 1 < C(l + m 2 /- 1 + mm (1+2 ™ ) ' 3 ' 1 / ^(r)y 1+ 2^ (r)dr) 

■/ o 



and we have: 

n-i— 2^i_ g i 1 
mm 2m = m^ -1 2m ^~ 1 '> . 

We have then: 

?r(i) < C(l + m^- 1 )^- 1 ) + m 2 ™^- 1 ) / ^(r)y 1+ 2^(r)«ir) 

JO 

where - ^ = (^-i^-i) < °- 

Recalling that (5 > 1 and </?(i) G L 1 (0,T) we find that for m big enough: 

V p (t)<C{C 1 + /%(r)/(r)dr) 

JO 

whence by Gronwall inequality: 

y(t)<C, tG[0,T], 
14 



where C depends on t. We thus have: 

\\p\\ L 2 m/ 3 + i < CmT^, t E [0,T]. 

Hence the inequality: 

\\p\\i*p)(t) <CkT&, t€[0,T\. (5.38) 

is valid for every k > 1, with C a positive constant independent of k > 1 but depending 
of the time. □ 

Remark 6 Let us point out that the estimate \5. 38\) is the key point in order to improve 
the range on f3. Indeed this last one is a refinement of the corresponding one in f32f . In 
particular we will be able to obtain the energy estimates \5. 66\) only with assuming (3 > 2. 

5.2 Second a priori estimate for the velocity 

In this section, we are going to furnish estimates on the velocity by using the gain of 
integrability on the density proved in the previous section. We are going essentially to 
follow the proof of [32j and to emphasize on the key point where we will only need the 
hypothesis (3 > 2. We begin with recalling some equation on the effective pressure defined 
in [23] and the rotational curl. We set: 

A = cmlu and B = (2 + X(p))divu - P(p), 
L = ~{d y A + 8 X B) and H =-{-d x A + d y B). 

We now want to obtain some estimates on the unknowns A and B, let us start with 
rewriting the momentum equation under the following eulerian form: 

dtu + u-Vu- -Au - -V((u + A(o))divw) + V(^^) = (5.39) 
P P IP 

Next if we apply the operator curl, we get: 

d t A + u ■ VA + Adivu = d y L - d x H. (5.40) 
Next we apply the operator div to the momentum equation f|5.39|) : 

d t divu + u-Vu- -Au - -V((p + A(p))divu) + V(^^) = 0, (5.41) 
P P IP 

and via the mass equation we have: 

80 + U.VB- P (2 + A)(B(^)' + (^ 1 )>v„ p 
+ (2 + A)(t/J + 2C/„Vi + If) = (2 + \){L X + H y ). 
As in [32j multiplying the equation (|5.40p by A and integrate over Q we obtain: 

/ ~[A 2 ]dx + \ [ divuA 2 dx+ [ (Ld y A- Hd x A)dx = 0. (5.43) 
Jn2dt 2 J n Jq 
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Similarly multiplying the equation (|5.42p by 2+\^ an( ^ integrate over O we have: 

[ 1 d ,1 If B 2 If . 1 , o , 

/ -— -B 2 <ix-- / divu -dx-- / u-V - — -)B 2 dx 

J a 2 + Xdt K 2 ' 2j n 2 + A 2j n V 2 + A ; 

- p Bdivu{B(^)' + (^)')^ + jf B(U 2 X + 2U y V x + V y 2 )dx (5.44) 

+ / (Ld x B + Hd y B)dx = 0. 
J n 



We recall now that: 



By combining the previous equality and (|5.44p we get: 

- [ P Bdwu(B(—^—)' + (^—)')dx + [ B(divu) 2 dx (5.45) 
Jn 2 + X 2 + A J n 

+ 2 f B{d y Ud x V - d x Ud y V)dx + f (Ld x B + Hd y B)dx = 0. 
Jn Jn 

Summing (|5.43[) and (|5.45p we have: 

Md r/lJ £ 2 , f 1 , „ 2 , /■ (A. + B^) 2 + (-A x + £„) 2 , 

- 1 y 5 2 dhm(^ - P(^)') + 2 | B(tr y ^ - E^ds 

- f pBdivu(B(- 1 —)' + (——)')dx + [ Bdivu 2 dx = 0. 
Jn 2 + X 2 + A Jq 



As: 



divu 2 = divu(^ + ^L), 



we deduce: 

IA^ 2l ^B 2 , fU, , 2 ._, A (A y + B x ) 2 + (-A x + B^ 



—dx 

n P 



IM [{A2 + ^ )]dX + L> UA2dX + In 

2 [ B{U y V x -U x V y )dx- [ pBdwu(B(- 1 —)' + (-^—)')dx = 0. 
Jn Jn i + a i + a 



(5.46) 



+ 

Let us set: 



, (Ay + B X ) 2 + (~A X + By) 2 , , I r 

a(t) = ( K y XJ v 2/7 dx) 2 , tE[0,T]. 



(5.47) 
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Next we have: 



( {{A y + B x f + {-A x + B y ) 2 )dx= I {A 2 x + A 2 y + B 2 x + B 2 y )dx. 
Jn Jn 

Let us observe that for every r, 1 > 4r > 0, from the result on elliptic system and by 
Holder inequalities we get as in [32J: 

\\VA\\ T , 2 ^ r un, + \\VB\\ r , (1 - r)f n,<C( ( ( A y + B *) 2 + i- A * + B yf dx) h 



iL^-^iQ) ' II v - D llL 2 ( 1 -'')(n) 

Jn r 

( / p^dx)^^. 

From (|5.38p . we have: 



n 



|VA|| L2(1 - r){n) + ||VS||^ (1 -r) ( n)) < C{-)^=Da{t) (5.48) 



Remark 7 Let us point out that the estimate fl5.^ff| ) is better than the corresponding one 
in due to the better estimate 115. 38]) . 

Moreover via (|5.38|) we also obtain the following inequality: 

(l|V«|| L2( n) + \\A\\ L 2 (n) + ||V2 + A divu\\ L 2 (n) ) < C(l + Z(t)), t e [0, T). (5.49) 
Now, are interested in providing other estimates for the non positive terms of the equality 

Estimates for the terms of (|5.46j) 

Following |32| . using (??), the lemma [T] (with a = 2 (i-2e) ) anc ^ Young's inequality (with 
P = 2 ^\Z 2 e £ \ Q = ^T=§r and pi = [izMg±g) ; Ql = the , we obtain: 

\\ ! div«^ 2 |<i||divn(t)|| L2 p(t)||i 4 , 
j n 

1 l-3e 1 1-e 

< C||divu|| L 2(2 ||t1|| L 2)i-2e (2 ||V^|| i 2(l- £ )) 1 - 2e 

< C(l + Z(t))Z(t)^({^)^)a(t))^ 

< 5a 2 (t) + C(<5)(1 + Z^)) 2 ^ 1 ^) 2 ^)/^ 

< <5a 2 (t) + C(5)(l + Z(t) 2 ) 2+ ^(i)^. 
We now are interested in estimating the term in (|5.46p corresponding to: 

J ' = li/„ B2di ™(2TA-"(2TA»>l 



1 ( n2 , B P w 1 , 1 ,,v, I 
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Easily there exist a positive constant C > such that: 



1 2 + A rv 2 + A' 
for all p 6 [0, +oo). We deduce that: 

By Young's inequality we have: 

Now, the last term in (|5.46p can be treated as follows: 

\2 [ B(U y V x -U x V y )dx\< [ \B\(U 2 + U 2 + V 2 + V 2 )dx 
Jn Jn 

Via the previous estimate, the notations (|5.47p . and using the equality (|5.46p . we get: 
i(Z 2 (t)) + a 2 (t) < 5a 2 (t) + C(5)(l + Z(t) 2 ) 1+ T^(-) T=i 

r |d|3 r (5.50) 

+ C ( 1 + 2 + A^ ) + / \B\(U 2 + U 2 + V 2 + V 2 )dx. 

It remains to estimate the two last terms on the right hand side of (|5.50p . In this goal, 
from (|2.1U|) we have: 

\\B\\ L ^ m < C(||fl|Ui (n) W||VB||^ ||S||^ 2(1 _ £)(n) ) (5.51) 

fl— e) 2 

where: s = _ n _ * and C > is a positive constant independent of m > 2. 

Now in inequalities (|5.50p and (|5.5ip we fix e = 2~ m with m > 2. Using estimate (|5.38|) 

for the density, we derive the inequalities: 

\\B\\ L i(n) = I \B\dx= I (^)l|B|(2 + A)lcZx<||(2 + A)l|| x2(n) Z(t) 
< CZ(t). 



( / (^-) 1 - £ | J B| 2 ( 1 - £ )(2 + A) 1 - £ dx)^) 
iji 2 + A 



IZJ|| ^-)(n) - vy^aTA 

< Z(t) s \\2 + Ap!_ e < C(-)^)Z s (t) < C2 2sm Z s (t), 
L—(fi) e 

< CZ s {t). 
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(5.52) 



From inequalities (|5.48p and (|5.5ip we finally obtain: 

\\B\\ L2m(n] < C(Z(t) + mk™)^(a(t)) 1 - s Z s (t)), 

< C(Z(t) + m^(j)^T)(a(t)) 1 - s Z s (t)). 

Now dealing with the integral with |i?| 3 , we have: 

I-BI 3 f \B\ 2 ~~^ 1 i !, i 

- L - L -drc= / 11 ; ( -)^^\B\ 1+ ~dx 

n 2 + A (2 + X) 1 2(m-i) 2 + A 

f \B\ m " 1 m n 1 -22- 

< / '' i \B\—^dx<Z(tf-—^\\B\\^ 

Jn (2 + A) 2(m-l) V ' 

^ {f^)^dx) l -^{fjB\ 2m dx)^ <~ 1+ W^ Z(t) 2 -^\\B\\^ (n) 

o 1 m m 771 m(l-s) ro(l-s) ms . 

< CZ 2 ~(t)(Z(t)~ + m 2 ^- 1 ) (a(t))~^^ Z~ (t)) , 

. „ m 777 m(l-a) m(l-s) „ ms-1 . 

< C(Z(tf + 777 2 (— !)( — )2(/3-l)(m- 1) (a(t))^^Z i+ ^T(t)), 

where C > is a positive constant independent of m > 2 and e = 2~ m . Finally applying 



applying Young's inequality with p = ^(i-s) and q = , m 2 ^ 1 ^2 2 we have: 



2 + A ~ v v ' v e 

n . „ m-l 777 m-l , 2(l-ms) 

< 5a 2 (t) + C(5)(Z 3 (t) + m ™( S +i)-2(_ )(/3-i)(m( s +i)-2) 2- 4+ ™( s +i)-2(t)) 
From (|5.5ip . we verify that: 



(5.53) 



m(l -e) 2 m(2-e) + e-l 

1 - 777S = 1 = e ^ '— _ 5.54 

m — e(l — e) m — e(l — e) 



and then: 

lim (2 m (l -ms)) = 2, 

rn— >+oo 

hence via (|5.53|) and (|5.54j) we get: 

l^l 3 ^/X„27^ , ^U,-, , ~2,^2 , „,"W 



< 5a 2 (t) + C(S)((1 + Z 2 {t)) 2 + m(— )^(1 + Z 2 (t)) 2+ ^+^) . (5.55) 



1 — ms 



2 + A ~ w wvv 
Now, consider the last term in (|5.50p : 



I 2 = I \B\(U 2 + U 2 + V x + V 2 )<\\B\\ L2 m {n) ([ (\VU\ 2 + \VV\ 2 )^dx) 1 -^. 
Jn Jn 

Recalling the relation 3 > 2 m^-i > 2, m > 2, from the properties of elliptic system [?] we 
derive the inequality: 



( / (|VC/| 2 + | Wl 2 )^^) 1 -^ < C(||div7z|| 2 4m 
Jn l^^t 



+ PII 2 4m ) 
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Thus the previous inequality furnish the estimate: 

h < C\\B\\ L2m(n) (\\divu\\] ; 4™ + ||Af 4™ ). (5.56) 

£2m-I (Q) £2m-l (Q) 

Next we have as ^4 vanishes on the boundary of the domain f2, we applythe Gagliardo- 
Niremberg inequality: 



l-e 



:(0) 

< CZ 2 ~ m ^-^) (t)((-)W^a(t)) m ' 1 - 2£ ) < CZ 2 "" 1 ' 1 ^) (*)(a(*)) 

Since 5 = (2 + A)divu — P, estimate (|5.38p provides: 

B P 

||divn|| 2 4m =11- r + rll 2 4m , 

<c(||-^-|| 2 *™ +1). 

2 + A L2m=T(fi) 

We can now deal with the right-hand side of (|5.58p as follows: 

2m(2m — 3) 

B P o / f \B\ (™-l)(2m-l) 2m . ,_J_ 



(5.57) 



(5.58) 



+ 7^11! 4m _ < ( / ' ' , , \B\(^)(^-l) dx y 



"2 + A 2 + X n L^^(n) ~ y J Q 2 + A 

1 r |_g|2 2m-3 1 /" \B\ 2 2m-3 

L 2m fH)( / (m-l)(2m-l) 

-/H (2 + ^) m(2m-3) Jf7 + A j 



- 1 - 2m-3 — 5— r, 1 

1 I ryi + W o_ o — II DM m-1 / -7 AA ;^3T 



Thus, 

||divn|| 2 4m <C{l + {Z(t)f-^\\B\Y^ ). (5.59) 
Using estimates (|5.57|) and (|5.59p , from (|5.56|) we have: 

h < C\\B\\ L2m{n) (l + Z 2 -^) (t)( a (t))TTO + Z 2 -^(t)||B||^ (n) ). 
Using estimate (|5.52p for ||-E> Hj^m^), we finally get: 



J 2 < C((Z(i)) 3 -d- 2 -) (a^))^ 1 - 2 -) + + Z 3 (t) 



+ m |(^)2^ Z 2 + s -^H7i( t )( a ( t )) 1 - s +MT^y 

1 ?7i m(l-s) n 1 ros-1 m(l-s) 1 ffl l-s 1 

+ m5 (— ) 2(/3-D(m-i) Z 2+ m=r(i)(a(t)) «-i + m5 (_)2</»-i) Z s (t)(a(t)) 1_s ) . 
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Using Young's inequality, we treat the summand in (|5.60p as follows: 

C{Z{t)f-^^) ( a (£))T^§iy < Sa 2 (t) + (7(1 + Z 2 (t)) 2 , 
C(Z(t) + Z 3 (t)) < C(l + Z 2 (t)) 2 , 

Cm2(— ) 2(/3-i) Z + M 1 -^) (i)(a(i)) s+ -d-2 £ ) < 

try) 1 q | 1 — ms+ (2ms — l)e 

5a 2 (t) + Cm(-)Fi(l + £2^+ (1+i)m(1 _ 2e) _ 1+e> 
Cm^{™)^^Z s {t)(a{t)) l - s < 5a 2 (t) + Cm(™)^(l + Z 2 (t)), 

1 771 m(l — s) d i ms-1 m(l — s) 

m 2( — )2(/3-l)(m-l) Z^ + ^^(t)(a(t))^^" < 

o 777, 1 r> o I 1 — ms 

5a 2 (t) + Cm(j)7=i(l + Z 2 (t)) 2+ (i+*>™-2 

Here <5 is a small positive constant to be mentioned below. From inequality (|5.6tjp we 
derive that: 



h < $a 2 {t) + C(m( j)~(l + Z 2 (t)) Z+ ^+s)ra(l-2e)-l+e 



rj i 1 — ms 



(5.61) 



+ (1 + Z 2 (t)) 2 + m(-) — (1 + Z 2 (t)) + m(-) — (1 + Z 2 (t))^a+^) 
From (|5,55p and (j5.61|) . and inequality (|5.50p we have: 
\j t (Z 2 (t) + a 2 (t)) < 6a 2 (t) + C(5)(l + Z 2 (t)) 2+ ^ (i) A 
+ C(l + <5a 2 (t)) + CC(<5)(1 + Z 2 {t)) 2 + CC(5)(m(-)Fi(l + Z 2 (t)) 2+ 

(5.62) 



e 

777 i r) I 1 — ms-\-{2ms — l)e 

45a 2 (t) + C(((l + Z 2 (t)) 2 + m (—)J=i(l + Z 2 {t)) 2+ a+.) m (i-2 e )-i+. 



e 

+ m (^)^T(i + z 2 (t)) 2 + m( ™)^(1 + ^ 2 (0) 2+m(ls + T) - 2 )- 
Choose (5 > such that: 

55 + 5C = ~. 

Since s = — — r. — \ and e = 2 171 , m > 2, we have: 

1 — ms , 1 — ms + (2ms — l)e e 

< 4e, - r — — r < 4e and < 4e. 



m(s + l)-2~ ' (1 + s)m(l - 2s) - 1 +e ~ 1 - 3e 

Then by (|5.62p and the fact that Z 2 (t) G L 1 (0, T), we obtain the inequality with < 
f <f: 

— {I + Z 2 (t)) + a 2 (t) < m(™)^(l + Z 2 (t)) 2+4 % (5.63) 
From ()5.63p we have for < i < T: 

1 1 ^ ,fn^-^— 

+ Cme{—)P- 1 > 0. 



(1 + Z 2 (i)) 4£ (1 + Z 2 (f)) 4e v e 
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Remark 8 Let us point out that the last inequality is better than in [32] and allows us 
to assume only (3 > 2. 

Now, take N > 2 such that: 

1 - CNe{— + Z 2 (0)) 4e > i e = 2~ N . 

S — 

Here the fact that /3 > 2 allows to conclude and by this fact improve the results of [32j . 
We get finally that for < t < T: 

Z 2 (t) < 2 2JV ~ 2 (1 + Z 2 (0)) - 1, te[0,T]. (5.64) 

Now, from inequality (|5.63p we get moreover that: 

(-T 

a 2 {t)dt < C. (5.65) 

o 

Now by estimate (|5.38p for the density, there exists a positive constant C depending 
continuously on the data of the problem and such that: 



sup / ( (curfu) 2 + \ ((2 + X(p))divu - P(p)) 2 ) (t,x)dx < C, 
<t<TJn\ 2 + A(p) J 

sup / ((curfu) 2 + (2 + A( / o))(divn) 2 ) (t, x)dx < C, (5.66) 
<t<T Jn 

dxdy < C. 



o<t<T Jn 
' T '■ (A y + B x ) 2 + (-A x + B y ) 2 



'o Jn P 

The rest of the proof follows exactly the same lines than in [32] and then we refer to [32] . 

6 Proof of theorem 11.21 

The existence part of the theorem is proved by an iterative method. We define a sequence 
(q n ,m n ) such that: 

' d t q° - pAq° + divm = 0, 
< dtm° — /jAm° + rm° = 0, 
k {q°,m°) = (q ,m ). 
Assuming that (q n ,m n ) is in Ej- with: 

E T = [QriBfr) n 4(^1 )) x (Ct(B 2 % ) n L x T {Bj+ n Bft )f , 

we define then g n+1 = g° + <f + , m n+1 = m° + fh n+l such that (g n +i, m n+ i) be the 
solution of the following system: 

' d t g n+1 - / uAg n+1 + divm n+1 = 0, 
< <9 t m n+1 - pAm n+1 + rrh n+1 = G n ; 
i (r +1 ,m"+ 1 ) A=0 = (0,0), 

with: 

771 71 

G n = -div(— ®m n ) 

We also set: h n = q n + 1. 
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1) First step, uniform bounds: 

Let e be a small parameter and by proposition 13.41 we have for any T > 0: 
WQ II M- M-+2 

L T (• B 2,l) nL T( B 2,l ) B 2,l ,„ „_ N 

n (6.67) 

llm II n , n _, n,, < C mn at,. 

L»(B^ )DLJ,(BT nB^ + ) b£ 
We are going to show by induction that for e > small enough: 

(p„) ii(r,m n )iiF T <e. 

As (foi^o) = (0)0) the result is true for n = 0. We suppose now (P n ) true and we are 
going to show (V n +l)- 

To begin with we are going to show that 1 + q n is positive. Indeed we have: h = h\ + h 2 
such that: 

d t h1 - nAhl = 0, 

and: 

d t h° 2 - ^iAh° 2 = -divm , 
(hi)/t=o = h . 

By proposition fj3.4f) and ()6.67p we have for any T > 0: 

\\h° 2 \\ N < C\\mo\\ n (6.68) 

By maximum principle, we have for any t > 0: 

h\(t,x) > min ^o(^) > c > 0. 

We deduce that for 77 = ||mo|| w_ 1 (at least inferior to ^ with the C of (|6.68p ) small 

B 2,l 

enough and any t > 0: 

h°(t,x)/geqj > 0, 



and 



Z°(t,x)>|-1. 



and by definition of g n and the assumption DP n that: 

<f(i,x)>^-l-e. 
In particular for e small enough at least e < § , we deduce that: 



/i" = 1 + g " > £ > 0. (6.69) 



23 



In order to bound (q n , m n ) in Et, we shall use proposition ^. 4l and in particular estimating 
G n in L^(i? 2 2 1 ). By using proposition 13.21 (|6.69p and lemmadl we obtain: 



VPtP 1 TVtP 1 

||div(— ®m n )\\ n , < \\— ®m n \\ n , , 
<C||m n || 2 ,v (||— ^ 1|| iv +1), 

<C(||m°|| 2 n +\\m n \\ N) 2 {\\q n \\ n + \\q°\\ N+l), 

Ll(B^) L^(B^) L^(B^) ' 



(6.70) 



Therefore by using (|6.70p . the proposition 13.41 and (V n ) we obtain for any T > 0: 

) (6.71) 

By choosing i] = e and e < 2C(2+||q C) jj * ms i m P nes C^On+i- We have shown by 

induction that (q n ,m n ) is uniformly bounded in Ft for any T > 0. 



^1 n+1)[ | < C (|| m 0||2 ^ +e )2( £+ || 0|| +1 ) 

< C(r ] + e) 2 (2+ \\q°" 



Second Step: Convergence of the sequence 

We shall prove that (q n ,m n ) is a Cauchy sequence in the Banach space Ft, hence con- 
verges to some (q,m) G Ft- 
Let: 

5q n = q n+l - q n and 6m n = m n+1 - m n . 
The system verified by (5q n ,5m n ) reads: 

' d t 5q n - fiA5q n + div<5m n = 0, 
d t 5m n - nA5m n + r5m n = G n -Gn-l, 
<5g n (0) = 0, <5u n (0) = 0. 

Applying propositions 13.41 and using (V n ), we get for any T > 0: 

\\(5q\5m n )\\ FT < C||G--G n _ 1 ||^ (BW/2 - 1) , 

By using proposition 13.21 and lemma [U we get: 

||(<Wm")|| FT < CelK^™- 1 ,^- 1 )!!^. 

So by taking e enough small we have proved that (q n ,m n ) is a Cauchy sequence in Ft 
which is a Banach space. It implies that (q n ,m n ) converge to a limit (q,m) in Ft- It is 
easy to verify that (q,m) is a solution of the system (|1.5p . 
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3)Uniqueness of the solution: 



The proof is similar to the proof of contraction, indeed we need the same type of estimates. 
Let us consider two solutions in Er- (gi,mi) and ^2,1^2) of the system (jl.5p with the 
same initial data. With no loss of generality, one can assume that (qi,mi) is the solution 
found in the previous section. We thus have: 

We note: 

Sq = q 2 — qi, dm = m 2 - mi, 

which verifies the system: 

dtSq — fiA5q + divJm = 0, 

dtSm — uA5m + r5m = — div( — — ® mi) + div(— — (g) mi) 

h\ hi 

By using proposition 13.21 13.41 and lemma [Hon [0, T±] with < T we have: 

\\(5q,5m)\\ ET < A(T)\\(5q,5m)\\ ET , 

such that for T small enough A(T) < |. We thus obtain: 5q = 0, Sm = on [0, T]. And 
we repeat the argument in order to prove that: 5q = 0, Sm = on M + . This conclude 
the proof of theorem 11.11 □ 
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